OKCon2010

r

Clear Climate Code

clearclimatecode.org

Nick Barnes
David Jones

Ravenbrook Limited

... and other contributors.

v

Ravenbrook

J

0.5

0.0

GISTEMP
Global Temperature Anomaly

0.5

-0.5

0.0

-0.5

1880 1890 1900 1910

10 1950 1960 1970 1980 1990 2000 2010

Ravenbrook

0.5

1910 1920 1930 1940 1950 1960 1970 1980 1990 2000 2010

Ravenbrook

J

r

FORTRAN 77

DO 20 JC=1,JCM

LT100S(JC)=NINT(100./PI180*ASIN(SNS+(JC-1)*DSLATJ))
LT100N(JC)=NINT(100./PI180*ASIN(SNS+JC*DSLATJ))
SNLATJ (JC)=SNS+(JC- .5)*DSLATJ
20 CSLATJ ( J C)=SQRT(1. -SNLATJ (JC)**2)
DO 50 I=1,NUMJ(J)
NR=NR+1

C**** Complete definition of SKIP-array
DO 30 IN=INFRST(J) , INLAST(J)
30 SKIP(IN,NR)=. FALSE.
DLON=360./NUMJ (J)
XEAST=-180.+I*DLON
XWEST=XEAST-DLON

C**** Extend each box by half a box in each direction (>Rcrit km ?)
XW(NR)=XWEST- .5*DL0N
IF (J .EQ.1.0R. J .EQ.8) XW(NR)=-180.
XE(NR)=XEAST+.5*DL0N
IF (J .EQ.1.0R.J .EQ.8) XE(NR)=180.
XN(NR)=ASIN(SNN)/PI180+DDLAT
XS(NR)=ASIN(SNS)/PI180-DDLAT

V

Ravenbrook

J

From the Internet

After a few weeks of effort no one was able to get the code to compile and execute.

Nobody has the patience to slog through the crap.

I no longer acknowlege GISTEMP as a legitimate source of tem[p]erature data

It is quite evident that NOAA, GISS, GISTEMP and GHCN are skewing their datasets toward the higher temperatures

One has to wonder though, just how much of that has to do with temperature measurements being taken in
proximity to a growing human population at these locations?

One would have a h[a]rd time discriminating between the novel bugs that did not affect the
published results from the prior bugs that did.

nobody (that I know of) has been able to get the FORTRAN code released by GISS last year to run.

Following that, an analysis of the methodology and replication of the computer program
output to see if it matches the current data sets

Grilling the Data
everybody knows adjustments are stupid
Gisstemp was like divinci code written by the shortbus crew

Ravenbrook

Motivation

ARE YOU COWNG TO BED?

J I CAN'T THIS
/ 15 IMPORTANT

WHAT? /

J SOMEONE. 15 WRONG
ON THE INTERNET.

I

r

Clear Climate Code Goals

1. To produce clear climate science software;

2. To encourage the production of clear climate science software;

3. To increase public confidence in climate science results.

v

Ravenbrook

J

Plan for GISTEMP

0. Run GISTEMP as it stands, step by step. Retain all intermediate
data files.

I. Reimplement each step in Python, with the same inputs and outputs,
using GISTEMP intermediate data files as test data.

>

2. Combine the Python steps into a single program, simplifying by
eliminating the need for intermediate data files.

3. Consider the correctness of the underlying algorithm.

4.Webify, data visualisation,

V

Ravenbrook

J

0.5

0.0

GISTEMP
Global Temperature Anomaly

0.5

-0.5

(and we match this)

0.0

-0.5

1880 1890 1900 1910

10 1950 1960 1970 1980 1990 2000 2010

Ravenbrook

Code Met

8461

misc

IliTtTi

3562

Python

GISS

ccc

rics

3562

comment
docstring

blank
code

gistemp

Ravenbrook

J

Independent Analyses

Global Temperature Reconstructions

NASA's response

"I hope to switch to your
version of that program ...
Ideally, we would like to replace
our whole code."

— Reto Ruedy, 2010-01-27

r

SciTech Committee

54. It is not standard practice in climate science and many other
fields to publish the raw data and the computer code in academic
papers. We think that this is problematic because climate science is a
matter of global importance and of public interest, and therefore the
quality and transparency of the science should be irreproachable.

... other countries have followed suit impeding the flow
of processed and raw data to and between researchers.
This is unfortunate and seems inconsistent with policies
of open access to data promoted elsewhere in
government.

v

Ravenbrook

J

Future

Clear

Your Climate Science Code?

imate Code Projects

0.6
0.4

Reconstructed Temperature

Medieval
Warm

2004*

Little Ice Age

1000 1100 1200 1300 1400 1500 1600 1700 1800 1900 2000
Global average temperature 1850-2007

, .„„. Based on Brohan et al 2006

S

rrj
ir

a

0.6
0 4
0.2
0.0
-0.2
-0.4
-0.6
-0.8

Annual average and 95% confidence range
Provisional annual average (200S) and 95% confidence range
Smoothed annual average and 95% confidence range

1850

Met Office HsJIey Centre

1900 1950

Source: wvw.metoffice.gov. uk/hadobs

2000

Crown Copyright 2008

Ravenbrook

ClearClimateCode.o

join us!

ravenbrook.com

(and others)
hire us!